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The induction of dc electronic transport in rigid and flexible trans-polyacetylene oligomers ac- 
cording to the tj vs. 2uj coherent control scenario is investigated using a quantum-classical mean 
field approximation. The approach involves running a large ensemble of mixed quantum-classical 
trajectories under the influence oi u! + 2u> laser flelds, and choosing the initial conditions by sampling 
the ground-state Wigner distribution function for the nuclei. The vibronic couplings are shown to 
change the mean single-particle spectrum, introduce ultrafast decoherence, and enhance intramolec- 
ular vibrational and electronic relaxation. Nevertheless, even in the presence of signiflcant couplings, 
limited coherent control of the electronic dynamics is still viable, the most promising route involv- 
ing the use of fs pulses with a duration that is comparable to the electronic dephasing time. The 
simulations offer a realistic description of the behavior of a simple coherent control scenario in a 
complex system, and provide a detailed account of the femtosecond photoinduced vibronic dynamics 
of a conjugated polymer. 



I. INTRODUCTION 

Conjugated polymers are of interest for their broad 
technological applications [l|, d, Q and because they are 
model systems that offer insight into the properties of 
soft organic and biological matter. As almost every pho- 
tochemical, photophysical, spectroscopic and charge and 
energy transfer process in these materials involves dy- 
namics of photoexcited states, the possibility to manip- 
ulate, at a molecular level, the dynamical properties of 
these excitations by means of lasers with well defined 
phases may have profound technological implications. 

As a first step toward this goal, we are interested in ma- 
nipulating the dynamics of photoexcited electrons along 
the backbone of a conjugated polymer using the prin- 
ciples of coherent control 4]. Specifically, we want to 
induce net dipoles along the material without introduc- 
ing a bias voltage. For this we apply zero-mean laser 
pulses with frequency components to and 2uj. Such fields 
are known to induce phase-controllable dipoles or cur- 
rents in anharmonic symmetric systems [5| even when 
they have a zero-temporal mean, a phenomenon that is 
referred to as laser-induced symmetry breaking. This rec- 
tification effect first appears in the third order nonlinear 
response of the system to the radiation field @ . At this 
order the system mixes the frequencies and harmonics of 
the uj + 2uj field, generating a phase-controllable zero har- 
monic (DC) component in the response. The DC compo- 
nent in the photoinduced dipoles is typically of the form 
(fi) ^ e^e2ui cos{(j)2uj - 2<j)uj), where enuj and (j)nuj denote 
the amplitude and phase of the nuj component of the 
field. Hence, simply by varying the relative phase be- 
tween the two incident lasers one can exert control over 
the magnitude and sign of the symmetry breaking. 

Laser-induced symmetry breaking 0, Q has been 



demonstrated in a wide variety of systems ranging from 
atoms to solid state samples. Experimentally, it has 
been implemented for generating anisotropy in atomic 
photoionization jTi], photocurrents in quantum wells [8|, 
intrinsic semiconductors @ and metal surfaces [lCl |. 
as well as directed diffusion in symmetric optical lat- 
tices [ll|. Theoretically, it has been studied for gener- 
ating transport in dop ed (T3 | and bulk semiconductors 
through interband [l3j and intraband [13 ] excitations, in 
graphene and carbon nanosheets [lB| and in molecular 
wires [lEIiS, [3j among others. The setup is of interest 
since, with current laser technology, it can be employed to 
generate controlled transport on a femtosecond timescale. 

Of major concern when using lasers to generate elec- 
tronic transport is the influence of the lattice dynam- 
ics on the rectification. The lattice can induce ultra- 
fast dephasing [H, [l^] processes that have deleterious 
effects on the control. Only a few attempts to quan- 
tify this effect exist. In all of them, the explicit depen- 
dence of the dynamics on the nuclear degrees of freedom 
is eliminated, and the effect of the lattice has been mod- 
eled through phenomenological relaxation [0| , stochastic 
forces [3] or thermal baths [l8|. Although formulations 
of this type have a considerable domain of applicabil- 
ity, they are not appropriate here since the electronic 
dynamics and spectroscopic observables in conjugated 
polymers depend on the detailed dynamics of the molec- 
ular backbone [ll], [l^, ■ This coupling between, and 
mutual influence of, electronic and vibrational degrees 
of freedom gives rise to the very rich photophysics of 
solitons, polarons and breathers and constitutes an im- 
portant distinction between "soft" materials and rigid 
solids 1, 22, 24] . Furthermore, the constant exchange 
of energy between the electronic and nuclear degrees of 
freedom during photoexcitation keeps the vibrations out 
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of equilibrium, rendering the thermal description inap- 
propriate. 

Here we numerically investigate the effect of vibronic 
couplings on the applicability of laser-induced symmetry 
breaking in 7r-conjugated systems by explicitly following 
the dynamics of both electronic and vibrational degrees 
of freedom. The simulations presented below provide a 
realistic description of the behavior of a simple laser con- 
trol scenario in a complex system, and offer a detailed 
account of the femtosecond vibronic dynamics of a con- 
jugated polymer. 

As a minimal microscopic model for conjugated poly- 
mers we adopt the Su-Schrieffer-Heeger (SSH) Hamilto- 
nian for trans-polyacetylene (PA) [l|,[2J|. The SSH model 
treats the polymer chain in terms of a one-dimensional 
tight-binding model in which the tt electrons are cou- 
pled to distortions in the polymer backbone by electron- 
vibrational interactions. It neglects quasiparticles inter- 
actions, assuming that they are relatively weak due to 
screening. Despite of its simplicity, the SSH Hamilto- 
nian has been remarkably successful in reproducing the 
band structure and the dynamics of excitations in PA. 
With it, we follow numerically the highly nonlinear cou- 
pled dynamics of electronic and vibrational degrees of 
freedom in neutral PA chains during and after photoex- 
citation with uj -\- 2uj laser pulses of varying frequency, 
width and intensity. The simulations a re p erformed in a 
mean-field (Ehrenfest) approximation |25l . [26l . [27l . l28j in 
which the nuclei are treated classically and the electrons 
quantum mechanically. Mean-field dynamics is the sim- 
plest mixed quantum-classical method that allows trans- 
fer of energy between quantum and classical coordinates 
with a proper conservation of energy [2^ IH, [l^, HI], 
and where transitions between instantaneous eigenstates 
are allowed. Allowing for change in the occupation of 
the electronic levels is crucial because the laser is con- 
stantly inducing electronic transitions and, additionally, 
electronic levels can approach one another closely during 
the dynamics and may lead to nonadiabatic transitions 
between electronic states. It is also a tractable method, 
as shown below, for treating the vibronic dynamics of 
multiple electronic states in the presence of a laser. 

In order to incorporate the effects of lattice fluctuations 
on the dynamics we follow the evolution of an ensemble 
of quantum-classical trajectories. The initial conditions 
are obtained by using importance sampling for the nu- 
clear Wigner phase space distribution in the harmonic 
approximation. In this way, the resulting averaged dy- 
namical observables retain relevant correlations between 
the electronic and vibrational degrees of freedom and are 
subject to electron-vibrational induced decoherence and 
relaxation. 

This analysis begins by briefly describing the model 
Hamiltonian (Sec. Ill A[) and deriving the equations of 
motion for the vibronic dynamics (Sec. Ill Bp . Then we 
turn to the important problem of choosing initial condi- 
tions for the evolution (Sec. Ill Cp and we introduce the 
geometric and spectroscopic observables that are used 



to monitor the dynamics (Sec. IIID|) . Our main results 
are presented in Sec. IIIII Specifically, after discussing 
the properties of the initial state fSec. IIII A]) , we investi- 
gate the typical dynamics of the chain generated by short 
(10 fs) and long (300 fs) pulses, and estimate the typical 
electronic dephasing times in PA oligomers (Sec. IIII Bl) . 
Then, we turn our attention to the effects that the lat- 
tice dynamics has on laser-induced symmetry breaking 
fSec. lIII C|) . The main results are summarized in Sec. lIVI 

II. MODEL AND METHODS 

A. The SSH model 

The SSH Hamiltonian fH, ^ models the PA oli gomer 
as a one-dimensional tight-binding chain, each site repre- 
senting a CH unit. The Hamiltonian for an iV-membered 
chain has the form: 

-ffsSH = -ffjr + -ffir-ph + -ffph, (1) 

where 

N-l 

-io XI XI i'^i+l,sCn,s + ci^sCn+l.s). (2) 
n=l s=±l 

describes the hopping of tt electrons along the chain with- 
out spin flip characterized by the lowest-order hopping 
integral to- The operator cj^ ^ (cn,s) creates (annihi- 
lates) a fermion in site n with spin s and satisfies the 
usual fermionic anticommutation relations {cn,s, s'} — 
5n,m5s,s' ■ The TT-electron-ion interaction term is given by 

N-l 

n=l s=±l 

(3) 

where u„ is the displacement of the n-th site in the x di- 
rection from the perfectly periodic position x = na, with 
a as the lattice constant. The operator iJ7r-ph couples the 
electronic states to the molecular geometry and provides 
a first-order correction to the hopping integral with a as 
the coupling constant. Finally, the nuclear Hamiltonian 
is taken to be 

^ph = E^ + yEK+i (4) 

n— 1 n— 1 

where M is the mass of the CH group, Pn is the mo- 
mentum conjugate to Un and K is an effective spring 
constant. 

The SSH model is an effective [l^l empirical model 
for non-interacting quasiparticles in PA. The electron- 
vibrational coupling and the hopping integral can be 
viewed as parameters in which screening and other high- 
energy effects have already been taken into account. 
The effect of the residual interactions that cannot be 
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accounted for by a simple renormalization of the one- 
electron Hamiltonian are completely neglected. 

Throughout this work, we use the standard SSH pa- 
rameters for PA [H]: a = 4.1 eV/A, K = 21 eV/A^, 
to = 2.5 eV, M = 1349.14 eV fsVA^ and a = 1.22 A. 
Results using this set of parameters agrees qualitatively 
well with experimentally determined properties. 

B. Photoinduced dynamics 



TABLE I: Parameters and labels defining the femtosecond 
laser pulses used — exp{ — {t — Tc)^/T^){e^cos{ujt+(f)^) + 
e2ui cos{2Ljt + (f)2ui)). Here l2ui is the intensity of the 2tj com- 
ponent at maximum field strength. 



Label 


(fs) 


Tw (fs) 


£2^ (V A-^) 


^UJ 1 £2ul 


/2^ (W cm"2) 


fl 


900 


300 


8.70 X 10"^ 


2.82 


1.0 X lO'' 


f2 


900 


300 


4.00 X 10-2 


2.82 


2.1 X 10"' 


f3 


50 


10 


8.70 X 10-3 


2.82 


1.0 X lO'' 


f4 


50 


10 


4.00 X 10-2 


2.82 


2.1 X 10^" 



To simulate excitation of the chain by a femtosec- 
ond laser pulse, an interaction with an electric field 
in dipole approximation is added to the SSH Hamil- 
tonian. The electric field is taken into account to all 
orders and the dynamics is followed in the mean-field 
(Ehrenfest) approximation. Our approach resembles 
that of Johansson et al. [1^ IIO] and Streitwolf [si'l, 
but we dynamically propagate an ensemble of quantum- 
classical trajectories rather than a single trajectory that 
starts from the optimal geometry. Single-trajectory ap- 
proaches [21, 22, 23, 29, 30, 31, 32] offer valuable insights 
into the time scales involved in the vibronic dynamics, 
but are insufficient for our purposes because including 
the decoherence effects requires averaging over the nu- 
clear degrees of freedom. Further, during this study we 
have observed that the vibronic dynamics of the SSH 
chain tends to exhibit deterministic chaos and, hence, 
individual trajectories become inaccurate for sufficiently 
long integration times. By contrast, results for an en- 
semble of trajectories are meaningful even for long inte- 
gration times, provided the shadowing theorem |3J, [s^l 
holds. 

In the presence of a radiation field, the total Hamilto- 
nian of the system assumes the form 

Hs{t)^Hssii + HE{t), (5) 

where HE(t) describes the interaction of the chain with 
an external electric field E{t) in the dipole approximation 
and Coulomb gauge, 

HE{t)^-{^le+^^^)E{t). (6) 

Here = ^\e\J2n,s^rici^s'^„,s and = +\e\J2n^n are 
the electronic and ionic dipole moments, Xn = na + u„ is 
the monomer position operator, a the lattice parameter 
and — |e| the unit electronic charge. The Coulomb gauge 
is consistent with our open chain boundaries. As a field 
we use the femtosecond uj + 2uj Gaussian pulses detailed 
in Table [T] at a variety of frequencies noted in the text. 

In the mean-field approximation the nuclei move classi- 
cally on a mean- field potential energy surface with forces 
given by [H, [M 113, S 

OUn 

where denote the antisymmetrized TV-electron 

wave function. The electrons are assumed to respond 



instantaneously to the nuclear motion and, hence, 
satisfies the time-dependent Schrodinger equation 

-i/elcc(M(i),t)|*W), (8) 

where iJcioc contains the electronic and the electron- 
nuclei interaction terms of the total Hamiltonian, and 
u — (mi,U2,-'' jUpf). Equations ([7]) and ([5]) define 
the Ehrenfest method. Feedback between the fast and 
slow degrees of freedom is incorporated in both direc- 
tions in an average self-consistent way. Note that the 
mean-field approximation avoids the expansion of the 
electronic wave function in terms of adiabatic basis func- 
tions. Equations ([7]) and ([5]) can therefore be integrated 
directly, making the implementation of the method par- 
ticularly simple. 

The electronic part of the total system Hamiltonian 
[Eq. ©] reads 

N-l 

Hclcc{t) = ^ [-to + a{Un+l - Un)](cJ,+ i_5C„,s 

n—l,s , ^ 

(9) 

^~^n,s^?^-t-l,s) 1^1 ^ ^ '^n^'n^s^n,s-E (t^ . 
n— l.s 

Since Hdccit) is a single particle operator, the electronic 
properties of the system are completely characterized by 
the electronic reduced density matrix, defined by 

Pn,m{t) = ^(^WIcJ^^.C^.^IvE-W). (10) 

s 

From Eq. ([8]) it follows that the dynamics of Pn.m{t) is 
governed by 

ift^Pn,m(t) = ^(*(t)|[ct,,C„^„i/elcc(i)]I^W) 

s 

" ^ ^ {_hm,m' (j^^ P 71, 7n' ^m',n(^)Pm'.7n(^)) ; 

(11) 

where hn,m{t) — {n, s\Heiecit)\m, s) are the single- 
particle matrix elements of Heiec{t), and \n,s) = c]^ jjO) 
where [0) is the vacuum state. In order to integrate this 
equation it is useful to employ an orbital decomposition 
for Pn,m{t), as described below. 
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Let |e, s) be the eigenorbitals of the system at prepa- 
ration time, defined by the eigenvalue relation ifeiec(i = 
0)|e, s) — e|e,s). Using this basis, the initial electronic 
reduced density matrix can be expressed as 

JV 

Pri,m(0) = V y'(e,s|n,s}(m,s|e',s} 

i^=is (12) 

(*(0)|<,c,,,,|vI/(0)), 

where c| ^ creates a fermion with spin s in the or- 
bital with energy e (i.e. |e,s) = ct_^|0)), and 
(^'(0)|c| gCe/^s|5'(0)} characterizes the initial electronic 
distribution among the single particle states. In writ- 
ing Eq. (|12p we have employed the basis transforma- 
tion function Cn,s — and its hermi- 
tian conjugate, which takes into account the fact that 
(e, s|n, s') — (e, s|n, s)ds,s' for the Hamiltonian under con- 
sideration. Upon time evolution we assume that Pn,m{t) 
maintains the form in Eq. (H^. That is. 



N 



= Y{e{t),s\n,s){m,s\e'{t),s)x 

€,e' = l s 

{^io)\cic,,^.s\^m. 



(13) 



The utility of this ansatz is that if the time-dependent 
orbitals \e{t),s) satisfy the single-particle Schrodinger 
equation 



ih—\e{t),s} ^ H,i,,{t)\e{t),s), 
at 



(14) 



with initial conditions \e(t = 0),s) = |e, s), the reduced 
density matrix automatically satisfies the correct equa- 
tion of motion [Eq. pT|) ]. 

With the exception of the simulations presented in 
Sec. IIII B 3| the initial electronic state |5'(0)) is taken 
to be a single Slater determinant for which 



{^{0)\cl,c,>^s\^{0)) = S,,,,f{e,s), 



(15) 



where /(e, s) is the initial distribution function that takes 
values or 1 depending on the initial occupation of each 
level with energy e and spin s. In this case, Pn,m{t) as- 
sumes the simplified form 



N 



Pn 



= im^^W' ^)^'"' sHt),s)f{e, s). (16) 



Henceforth we drop spin labels since orbitals with oppo- 
site spin satisfy the identical equation of motion. 

Given Eqs. (O and pO]) and the total Hamiltonian of 
the system Eq. ([5]), the equations for the nuclear trajec- 
tories are given as 



Unit) 
Pn{t) 



_Pn{t)_. 
M ' 

-- - K{2Un{t) - Un+l{t) - Un-l{t)) (17) 

+ 2aRe{pn,„+i{t) -p„,„-i(t)} 
- \e\Eit) (p„,„(i) - 1) . 



The chain is taken to be clamped so that ui{t) = ujv(i) = 
andpi(t) = pN{t) = for all time, and Eq. (fT7|) is valid 
for n = 2,---,iV — 1. In turn, the orbitals that form 
Pnm{t) satisfy Eq. (fH)) . so that 

ih^{n\e{t)) = [-ta + a(u„+i(i) - w„(t))] (n + l\e{t)) 
at 

+ [-to + a{un{t) - Un-i{t))] {n - l\e{t)) 
+ \e\E{t) {na + Un{t)){n\e{t)), 

(18) 

for n,e = I,-- - ,N. Since the electrons are confined 
within the chain, {n\e{t)) — ior n ^ {1, • • • , N}. Equa- 
tions ini) and ([HI) constitute a closed set of N{N + 2) 
coupled first-order differential equations, which are inte- 
grated using the eighth-order Runge-Kutta method with 
step-size control [35.] . 



C. The initial conditions 

A crucial aspect of the problem is the choice of initial 
conditions. The oligomer is assumed to be in the ground 
electron-vibrational state. Hence, we first determine the 
optimal geometry of an A^-membered oligomer with open 
boundaries and clamped ends (ui = ~ 0) by an itera- 
tive self-consistent procedure, and then perform a normal 
mode analysis in the ground electronic surface as detailed 
in Ref. [H. This procedure provides the nuclear ground- 
state wave function in the harmonic approximation. A 
phase-space like description of the resulting nuclear quan- 
tum state is obtained by constructing the associated nu- 
clear Wigner phase-space distribution [33, [sl] function 
pw{u,p), which is just the product of the Wigner distri- 
butions associated with each vibrational mode 



N-2 



(19) 



Here Qj{u) is the normal mode coordinate of the j-th 
mode and Pj{p) its conjugate momentum. In the ground 
state, the Wigner distribution of each normal mode is 
given by [s^j 

PjiQj^Pj) = \eyiv{-'MLOjQyh)eyiY>{-Pf/nwjM), 

(20) 

where uuj is the frequency of the jih mode. The 2A^ — 4 



dimensional phase-space distribution in Eq. (|19p com- 
pletely characterizes the initial quantum state of the nu- 
clei. 

The ensemble of lattice initial conditions, 
{m*(0),p*(0)}, for the quantum-classical dynamics 
is obtained from a Monte Carlo sampling of the nuclear 
Wigner phase space distribution of Eq. p9|) . The aver- 
age classical energy of the resulting ensemble coincides 
numerically with the zero-point energy of the lattice. 
The associated initial values for the orbitals {(ri|e*(0))} 
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are obtained by diagonalizing the electronic part of the 
Hamiltonian in the initial lattice geometries {m*(0)}. 
Each initial condition i, together with the equations 
of motion pT|) and defines a quantum-classical 

trajectory (m*(0), p*(0), |vl^^(0)}) ^ (u^ {t) , (t) , [^^(t))). 
The set is propagated using a parallel algorithm and 
used to obtain ensemble averages. In this manner, 
the dynamics reflects the effects of lattice fluctuations 
and the initial quantum phase-space distribution of the 
nuclei. 

Since laser-induced symmetry breaking is an effect that 
depends on the third-order response of the system to the 
field, adequate convergence of the results with the num- 
ber of initial conditions requires a large number of initial 
configurations, 0(10^ x (2A^ — 4)). This is the main com- 
putational bottleneck of the present approach and limits 
our analysis to modestly sized oligomer chains. 

D. Dynamical observables 

During the vibronic dynamics induced by the laser field 
the molecular geometry distorts and this, in turn, induces 
strong changes in the electronic wave function. Below 
we describe the geometric and spectroscopic observables 
that we use to follow the evolution of both nuclei and 
electrons. 

Geometrical changes in the polymer backbone are 
characterized by the bond length alternation parameter 
r„, which quantifies the homogeneity in the distribution 
of TT electrons over the bonds. The ensemble average of 
this quantity is defined by 

'-^ (21) 

i—l 

where is the bond length between the n-th and {n + \)- 
th atoms along the chain in the «-th member of the en- 
semble and is the cardinality of the ensemble. When 
the alternation between single and double bonds is per- 
fect, the bond length alternation is constant (apart from 
end effects) and takes a value of tq w 0.08 A. An enhance- 
ment in the electronic delocalization tends to equalize the 
bond lengths in the polymer resulting in \rn\/r{) <C 1. 

The state of the electronic degrees of freedom is com- 
pletely characterized by the ensemble-averaged electronic 
reduced density matrix, 

1 ^ 

Pn,mit) = j;iY.Plrn{t) (22) 

i—l 

where „j is defined by Eq. pO]) . The diagonal elements 
of Pn m (i ) represent the density of charge along the chain. 
The off-diagonal elements represent the electronic coher- 
ences between different sites. 



A compact description of the electronic dynamics is 
offered by the polarization of the chain, defined by 

, , M N 
i—l n—1 

where xl^{t) = {na+u''^{t)) is the position of site n at time 
t in the i-th trajectory. The first term in Eq. P5)) comes 
from the dipole due to the nuclei, while the second one 
quantifies the electronic contributions. Any laser-induced 
symmetry breaking must manifest as a DC component in 
{fi{t)). In order to isolate the zero-frequency component 
from the oscillatory terms we integrate the signal over 
time 

iCit)}^-^ fdt'W)). (24) 

This gives the net cumulative dipole induced by the field, 
weighted over the chain length L and the temporal pulse 
width Tw. 

For each member of the ensemble the instantaneous 
eigenorbitals of the electronic Hamiltonian (including the 
field) are defined by the eigenvalue relation 

Hl,,Ml\t))-^,{t)\f{t)), (25) 

where Hl^^^{t) is the electronic Hamiltonian for trajec- 
tory i at time t. At each instant in time the instanta- 
neous eigenorbitals for different trajectories differ since 
each trajectory experiences a different vibronic evolu- 
tion. For each trajectory, at a given t, the orbitals 
e*(t)) satisfying Eq. (fT8|) may be expanded as a lin- 
ear combination of the instantaneous eigenfunctions I7'), 
k'W) = E^«;,.Wl7nO), where a\^,[t) = {iKt)\e\t)). 
The occupation number nl^ of the instantaneous states is 
thus given by 

(26) 

and varies as the system evolves. This contrasts with 
purely adiabatic dynamics with fixed level occupation. 
For the ensemble, we define the mean occupation and 
mean energy of the instantaneous eigenstates as 

(^''W) = A^E4(^)' K(i)) = ;^E<(*)' (27) 

i—l i—l 

for 7 = 1, • • • , iV. These quantities are physically accept- 
able if the ordering of the eigenstates does not vary 
among trajectories, a property that we assume. Through- 
out the text the label 7 is employed to number the in- 
stantaneous eigenstates of the oligomer in ascending or- 
der, with 7 = l,N/2 and N/2 + 1 denoting the lowest 
energy, HOMO and LUMO orbitals respectively. 
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III. RESULTS AND DISCUSSION 

Consider the ability ofuj+2uj fields to induce symmetry 
breaking in PA. The dynamics is characterized using the 
observables discussed in Sec. Ill Dl as well as by the differ- 
ent contributions to the chain energy. Vibrational effects 
are made explicit by contrasting the vibronic dynamics 
of the oligomers with the evolution of a single trajectory 
for an equivalent but rigid system initially in the optimal 
geometry. The chain is made rigid by multiplying the 
mass of the CH groups by a factor of 10®. In this way 
the lattice is made to move a thousand times more slowly 
while maintaining the same value of the electron-ion in- 
teraction. 



A. Initial state 

We study neutral oligomers with 20 sites and 20 tt 
electrons (L « 23 A) initially in the ground electronic 
configuration. This size was selected because it repro- 
duces qualitatively well the spectrum and dimerization 
pattern of longer chains without making the computa- 
tional effort prohibitive. The initial state for the flexible 
chain consists of an ensemble of 100,000 initial conditions 
obtained by sampling the ground-state nuclear Wigner 
distribution fSec. HTC)) . 

The resulting electronic and geometrical properties of 
the chain are shown in Fig. [TJ The initial geometry of the 
chain (Fig. \Vp) consists of a perfect alternation of dou- 
ble and single bonds, yielding (r„} ~ 0.08 A. The open- 
ends boundary condition results in a stronger dimer- 
ization near the edges of the oligomer. The structure 
is centro-symmetric with the inversion center residing 
between sites 10 and 11. The single-particle spectrum 
(Fig. [Tp and Table |TT| has a total width of Ato = 10 eV. 
It consists of N/2 fully occupied tt (valence) states and 
N/2, initially empty, tt* (conduction) states, separated 
by an energy gap 2 A = 1.8 eV. The reduced density ma- 
trix (Fig.[TJ'V-B) is concentrated along the diagonal of the 
plot reflecting the bonding pattern of the oligomer, with 
the spatial electronic coherence (measured by the magni- 
tude of the off-diagonal matrix elements) being ~ 10 — 15 
units. 

Note that the average geometrical and spectroscopic 
properties of the ensemble coincide with the ones ob- 
tained from a single configuration in the optimal geom- 
etry. However, within the ensemble there is a distribu- 
tion of values for the properties, displaying important 
deviations from the average. For example, we find that 
the ensemble of orbital energies form bands of states dis- 
tributed in a Gaussian fashion about the band averages. 
The banded-Gaussian distribution in energy space is a 
direct result of the Gaussian distribution of the initial 
geometries. The band-averages (e^) and bandwidths ai 
are enumerated in Table |TT1 Note that the bandwidth of 
the electronic levels near the energy gap (i = 10 — 11) 
is considerably broader than the one displayed by states 
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FIG. 1: Nuclear and electronic properties of a 20 atoms -|- 
20 TT electrons SSH-type chain in the ground-state configu- 
ration. Properties of the chain in the optimal geometry are 
compared with the average results obtained from an ensemble 
of 100,000 initial configurations. The upper panels show the 
reduced density matrix (chosen to be real) for (A) the optimal 
geometry and (B) the ensemble of initial states. The color bar 
is given in the far right. Panel (C): bond length alternation 
along the chain and (D) the single-particle spectrum, (o) - 
the optimal geometry; (x) - the ensemble average. In both 
cases, the initial average momentum for each atom is zero. 



TABLE II: Initial single particle spectrum for the 20 site neu- 
tral PA oligomer. The energies el were obtained at the opti- 
mal geometry, while (ei) denote the average orbital energies 
for the ensemble. The ensemble of orbital energies form en- 
ergy bands distributed in a Gaussian manner about the band 
centers, /(ti) = ^^^^ exp (— (e — (ei))^ /2af) , where at is the 
standard deviation. 
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-4.930 
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0.130 
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-4.735 


-4.755 


0.039 


12 
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-4.491 


0.046 
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2.045 
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0.091 
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0.107 
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4.735 


4.755 


0.039 


10 


-0.914 


-0.883 


0.130 


20 


4.893 


4.930 


0.034 



near the band edges. 

Last, due to residual anharnionicities around the equi- 
librium geometry, the distribution in Eq. (|19p is not com- 
pletely stationary. Rather, the average values of the ini- 
tial state display some fluctuations under free evolution. 
These fluctuations are very small, introducing only small 
changes to Fig. [H of the order of 1% in the band-gap 
energy and of 3% in the bond length alternation. 
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B. Rigid vs. flexible chain dynamics 

For rigid chains, changes in the electronic spectrum 
can only arise from the Stark shifts induced by the laser 
field. By contrast, for flexible chains there is a continuous 
exchange of energy between the electronic and nuclear 
degrees of freedom during the photoexcitation process. 
Below we present two illustrative examples of the pho- 
toinduced dynamics in PA. First we consider the case 
in which the lattice is photoexcited with a 10 fs pulse 
(Sec. IIII B ip . For such short pulses very little excited 
state dynamics occurs during the pulse. In the wake of 
the pulse we observe coherent vibrational breathing mo- 
tion that decays in ^ 200 fs due to vibrational energy 
redistribution. We then consider excitation with a 300 
fs pulse (Sec. HUB 21) . In this case the laser does not 
generate a breathing motion of the lattice. Rather, we 
observe a gradual increase in the electronic delocaliza- 
tion, accompanied by a concurrent change in the average 
spectroscopic and geometrical parameters of the chain. 

Last, in Sec. IIII iTS] we quantify the time-scale for elec- 
tronic dephasing for chains of different lengths. 
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FIG. 2: Energy contributions during the dynamics of PA 
under the influence of pulse f4 with hjjj — 1.3 eV and 
(I>2lj — 2(f)^ — for the rigid (gray lines) and flexible (black 
lines) chain. Panel (A): the electronic kinetic energy (includ- 
ing electron-ion interactions), (B): purely nuclear energy, and 
(C): radiation-matter interaction. The total energy is shown 
in (D). 

The different contributions to the rigid and flexible 
chain energy during and after photoexcitation with a 10 
fs (jj + 2ijj pulse (laser field f4 in Table|T| with hcu — 1.3 eV 
are shown in Fig. [2] Initially, the flexible chain has 0.83 
eV more energy than its rigid counterpart due to zero- 
point motion in the ground electronic state. Secondly, 
after photoexcitation, the rigid chain gains 1.85 eV from 
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FIG. 3: Occupation numbers of the instantaneous eigenstates 
near the band gap for the (A) rigid and (B) flexible chain 
during the dynamics induced by pulse f4 with hu = 1.3 eV 
and (j}2uj — '2(j>ui = 0. 



the laser while the average energy absorbed by the ensem- 
ble of flexible chains is 3.9 eV. This enhanced absorption 
of energy of the flexible chain is due to spectral broad- 
ening of states which form bands, with the broadening 
being especially large near the band gap (recall Table ITT)) . 
As exemplified in Fig. [31 this results in a wider set of elec- 
tronic levels being populated during photoexcitation. In 
both cases, the occupation numbers and the single par- 
ticle spectrum display the electron-hole symmetry of the 
SSH Hamiltonian. 

Figure d] shows the evolution of the bond length al- 
ternation parameter (upper panel) and the mean instan- 
taneous orbital energies for states near the energy gap 
(lower panel) of the flexible chain. The photoexcitation 
of the chain initiates a complex highly nonlinear vibronic 
evolution. The lattice is set to vibrate and large scale os- 
cillations in the bond length alternation are observed. 
By i = 67 fs the chain has reversed the sense of its 
dimerization and establishes a coherent "breathing" os- 
cillation of the bond length alternation in which (r„(i)) 
oscillates between the two senses of dimerization. Several 
carbon-carbon stretching vibrational modes with periods 
between 32 — 40 fs enter into the initial lattice dynam- 
ics. The breathing motion, however, is not long lived 
and the amplitude of the oscillations display an initial 
fast exponential decay with a characteristic time-scale of 

90 — 130 fs, followed by a slower damping. By t ~ 400 
fs most of the vibrational modes initially excited have 
already decayed, and gradually the system approaches a 
state of internal equilibrium through internal vibrational 
energy redistribution. 

These changes in the lattice geometry are reflected as 
changes in the single-particle spectrum. Upon initial pho- 
toexcitation the two states near the energy gap come very 
close to each other in less than 20 fs and the band-gap en- 
ergy fluctuates following the breathing oscillation of the 
bond-length alternation. A dip (crest) in the bond-length 
alternation corresponds to a minimum (maximum) in the 
gap energy. Concurrent with the dimerization pattern, 
the energy gap exhibits strong initial oscillations at time 
scales between 30-40 fs that rapidly decay, followed by a 
slower damping. The band-gap energy eventually reaches 
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FIG. 4: Vibronic dynamics of PA induced by pulse f4 with huj = 1.3 eV and (f>2oj — 20^; = 0. The upper panel shows contours 
of the bond length alternation parameter [defined in Eq. (|2ip ]. The lower panel, the mean instantaneous energies for the 
single-particle states near the energy gap. 
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FIG. 5: Decoherence during the dynamics of the flexible chain 
induced by pulse f4 (10 fs wide, centered at 50 fs) with fitj = 
1.3 eV and (j>2ui — 2(f)uj = 0. 



an equilibrium value of ~ 0.6 eV. The remaining energy 
levels also oscillate on time-scales determined by the dif- 
ferent vibrational modes of the chain. We observe that 
states around the band-gap couple to vibrational modes 
whose period is ^38 fs, while states near the band edges 
couple to lower frequencies modes. 

An estimate of the electronic decoherence that oc- 
curs during and after the pulse is provided by measur- 
ing the purity [H] Trp^ = J2„^rnPn,mPni.n, where p„ „, 
are the matrix elements of the ensemble-averaged elec- 
tronic reduced density matrix [Eq. (|22p ]. For pure sys- 
tems with a Pn,m{t) of the form in Eq. as is the 
case of the rigid chain, Trpp^^(, is a constant equal to 27V 
where Af is the number of electrons. For mixed states 
TrPmixod — Trppy^c and its decay offers insight into the 
decoherence time scales introduced by the coupling to the 



vibrational modes. Figure [5] shows the evolution of Trp^ 
for the flexible chain. The system begins in a state that 
is almost pure with Trp^/Trpp^^,, = 0.997. Upon pho- 
toexcitation the purity displays a sharp drop at ^ 40 fs, 
followed by some modulation in time scales commensu- 
rate with the vibrational periodicities of the chain. These 
modulations decay rapidly and the purity stabilizes at 
0.89 after ~ 400 fs. 

The dynamical behavior exemplified by our simula- 
tions describes quite well the main trends observed in 
recent sub-10 fs experiments that have been performed 
on PA [4^ and PPV [4lj samples: initial fast decay of the 
breathing motion (order of 200 fs) followed by a slower 
damping. Previous studies of the excited state dynamics 
in PA did not capture the decay since they only consider 
the evolution of a single trajectory. However, the time 
scales that we observe for the energy gap fluctuations 
(32-38 fs) are somewhat larger than the ones reported 
for PA of 31 fs and 23 fs. This is a consequence of the 
parameters employed for the Hamiltonian. 

We have observed that the basic features of the vi- 
bronic dynamics described above also apply to longer 
oligomers. However, in the long chain limit the two 
states near the energy gap can become degenerate, form- 
ing soliton-antisoliton pairs f23|, unlike to the case de- 
scribed above. 



2. The dynamics induced by a 300 fs pulse 

Consider now dynamics under the longer (300 fs) pulse 
fl with hu> = 1.18 eV. Here the pulse time width is 
long compared with the typical optical-phonon vibra- 
tional modes. Figure [S] and [7] show the different energy 
contributions to the chain during the dynamics and the 
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FIG. 6: Energy contributions during the dynamics of PA 
under the influence of pulse fl with hu = 1.18 eV and 
(j>2uj — '24'lo = 0. Panel (A) electronic kinetic energy (includ- 
ing electron-ion interactions), (B) the purely nuclear energy, 
and (C) the radiation-matter interaction. The total energy 
is shown in (D). The gray and black lines correspond to the 
rigid and flexible chains, respectively. 




FIG. 7: Occupation numbers of the instantaneous eigenstates 
near the energy gap during the dynamics induced by pulse fl 
with hu! = 1.3 eV and 02t^ — '24>ui = 0. (A) The rigid chain; 
(B) flexible case. 

instantaneous occupation numbers for states near the en- 
ergy gap, respectively. As before, due to broad nature of 
its transition frequencies the flexible chain absorbs 1.36 
eV more energy from the field, with a wider range of elec- 
tronic levels being populated during the process. The 
large-scale oscillations observed previously in the aver- 
age electronic and nuclear energy for the flexible chain 
(Fig. [^t^-B) are no longer present. 

The origin of this difference in behavior is readily iden- 
tified by considering the bond length alternation and 
the single particle spectrum during the dynamics (Fig- 
ure[8|) . During the first 400 fs the bond-length alternation 
and the transition frequencies display minor oscillations. 
These oscillations arise because our initial state is not 
an exact eigenstate of the Hamiltonian due to ground 
state residual lattice anharmonicities. Upon photoex- 
citation the bond lengths in the chain slowly equalize 
in order to improve electronic delocalization along the 
chain that permits accommodating the electrons that 



have been photoexcited. After the pulse the bond length 
alternation settles at ~ O.OSA. These changes in the lat- 
tice geometry are accompanied by a 0.7 eV red shift of 
the energy gap. However, contrary to the dynamics under 
the 10 fs pulse, the lattice does not establish a breath- 
ing motion. Instead the mean geometrical observables 
and spectroscopic quantities display a slow and steady 
change from the perfectly dimerized lattice to a chain 
with increased electronic delocalization. In essence, the 
frequency width of the pulse is not broad enough to bring 
about a coherent breathing motion of the chain. The 
latter requires participation of many vibrational states. 
These features of the vibronic dynamics induced by 300 
fs pulses also persist when considering longer oligomers. 

Decoherence occurs slowly and steadily during the time 
that the system is being driven by the laser field (Fig.lH]). 
No modulations of Trp'^ or further decoherence are ob- 
served after the pulse, and Trp^/Trpp^^,^ stabilizes at 
~ 0.88. Computations discussed below show that the 
electronic dephasing time is of the order of 2.5 fs. Hence, 
the limiting step in the decoherence process is not the 
rate at which lattice-induced decoherence occurs, but the 
rate at which the field is able to establish new superpo- 
sition states that then rapidly decohere. 

3. How fast is the electronic dephasing? 

A measure of the timescale for the electronic dephas- 
ing induced by intramolecular vibrational motions can be 
obtained by following the field-free evolution of neutral 
chains initially in an electronic superposition state. In a 
full quantum mechanical analysis, nuclear evolution on 
alternative electronic potential energy surfaces leads, in 
general, to electronic coherence loss. In our quantum- 
classical picture the electronic dephasing is captured by 
the distribution and evolution of orbital energies con- 
tained within the initial ensemble (recall Table |TT| . 

A detailed calculation reveals the features of this de- 
phasing. Figure [TO] shows the time-dependence of the 
polarization for PA oligomers initially prepared in a su- 
perposition l^'(O)) = ^(|G) + |E)), where |G) and [E) 
are the A^-particle ground and first-excited electronic 
states. The initial nuclear distribution is taken to be 
the one obtained in the ground electronic surface, and 
E) = clii^i g/'^e" ,s'\G) is generated by instantaneously 
promoting an electron from the HOMO to the LUMO 
orbital of the chain, so that e" ~ N/2. For chains 
with 4 sites the polarization displays a fast initial decay, 
and observes recurrences every ~ 30 fs. In a full quan- 
tum mechanical analysis these recurrences arise from the 
time dependence of the overlap of the nuclear wavefunc- 
tions in the ground and excited electronic states. In this 
quantum-classical picture, the recurrences are captured 
by the time-dependence of the orbital energies during the 
dynamics. Note that between consecutive recurrences the 
amplitude of the polarization gets reduced and {fJ-{t)) dies 
in hundreds of femtoseconds. For longer chains these re- 
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FIG. 8: The vibronic dynamics of PA induced by pulse fl with hu = 1.18 eV and (j)2uj — = 0. The upper panel shows 
contours of the bond length alternation parameter (defined in Eq. (|2H) ). The lower panel shows the mean instantaneous energies 
for the single-particle states near the energy gap. 
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FIG. 9: Decoherence during the dynamics of the flexible chain 
induced by pulse fl (300 fs wide, centered at 900 fs) with 
huj = 1.18 eV and (l)2u. - 2cj>^ = 0. 



currences are less important and the dephasing time is 
basically determined by the initial decay of the polariza- 
tion which occurs in less than 10 fs, with a characteristic 
timescale of 2.5 fs. These time-scales are in agreement 
with electronic dephasing timescales that have been de- 
termined in other systems using a fully quantum treat- 
ment El. 



C. Laser-induced symmetry breaking 

In the previous sections we identified important ways 
in which the vibronic couplings influence the photoin- 
duced dynamics of 7r-conjugated systems. The coupling 
introduces: (i) broadening of the electronic transitions 
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FIG. 10: Time dependence of the polarization {fJ.{t)) for neu- 
tral flexible PA chains with N sites. The initial state is a su- 
perposition, with equal coefficients, between the ground and 
first excited electronic state. The nuclei are taken to be ini- 
tially in the ground-state configuration. 



(see Table im and Figs. [3] and [7|), (ii) pronounced changes 
in the mean single-particle spectrum (see Figs. |4] and [8]), 
(iii) internal relaxation mechanisms, and (iv) ultrafast 
dephasing (see Figs. [5l iQland fTO)) . Because of this, when 
using lasers to control a dynamical process in this class of 
systems, the laser frequencies may become detuned from 
the desired transition and the control may be diminished 
by broadening of the electronic levels and internal relax- 
ations. Here we investigate the extent to which these 
additional complexities affect the ability to induce sym- 
metry breaking using the to vs. 2uj coherent control see- 
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FIG. 11: First moment of the charge distribution during and after photoexcitation with a symmetry breaking a; + 2lj laser 
field. (A) the flexible and (B) the rigid chain under the influence of pulse fl with hu — 1.18 eV and 021^ — 2(j)uj ~ 0. Panel (C) 
the resultant cumulative dipoles {C{t)) (black line, flexible chain; gray line, rigid case). Parts (A')-(C') as in (A)-(C) but for 
the f4 pulse with hcu = 1.3 eV and (l}2ui — 20i^ — 0. 
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FIG. 12: Frequency dependence of the degree of symmetry breaking generated by an + 2(jj field in rigid (gray) and flexible 
(black) PA oligomers with 20 sites. The initial energy gap of the oligomer is 2A = 1.8 eV. The plots show the asymptotic 
cumulative dipole (C(oo)) after photoexcitation with pulses (A) fl, (B) f2, (C) f3 and (D) f4, using (j)2ui — 2(j>uj = 0. 



nario. 

The results presented in this section focus on effects 
that depend on the relative phase of the laser field. This 
is because symmetry breaking contributions that depend 
on the carrier envelope phase are typically difficult to 
control (although not impossible, see e.g., Ref.lZsj) since it 
requires an experimental arrangement that permits lock- 
ing the carrier envelope phase of the incident radiation 
and that controls the position of the center of mass of the 
molecule with respect to the laboratory reference frame. 
This contrasts with relative-phase control that is unaf- 
fected by the center of mass motion and only requires 



manipulating the relative phase between the two central 
frequency components of the uj + 2u! field. 

We consider the evolution of rigid and flexible chains 
under the influence of the laser pulses in Table |T1 for 
different laser frequencies. The parameters chosen for the 
laser pulses encompass four illustrative cases: dynamics 
induced by weak and moderately strong pulses with time 
envelopes that are either short or long with respect to 
the typical 30-40 fs optical-phonon vibrational period. 
Specifically, pulses fl and f2 induce weak and strong- 
field dynamics, respectively, when the time-width of the 
pulse is relatively long (300 fs). By contrast, pulses f3 
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FIG. 13: The degree of symmetry breaking induced by an 
ulj + 2u field in rigid (gray) and flexible (points) PA, as a func- 
tion of the relative phase (j)2ui — 20;^ of the laser. The figures 
show the mean cumulative dipole {C(oo)) after photoexcita- 
tion with: (A) pulse f2 with ftoj — 0.9 eV; (B) pulse f4 with 
fujo = 1.3 eV; and (C) pulse f4 with hui = 1.6 eV. 



and f4 offer insight into the weak- and strong-field control 
when the laser pulses are short (10 fs), and of the order of 
the electronic dephasing time, permitting only a limited 
degree of excited state dynamics during the pulse. 

Consider first the dynamics in the weak- field/long- 
pulse (fl) and strong- field/short- pulse (f4) regimes. 
The relevant observable here is the polarization of the 
chain [Eq. which quantifies the symmetry 

breaking. Figure [11] shows representative results. Ini- 
tially (/i) — 0, reflecting the symmetric distribution of 
electrons over the chain. As soon as the laser is turned 
on, dipoles are induced in the system, oscillating with 
the frequencies of the field and various harmonics. The 
harmonic of interest is the zero-frequency (DC) compo- 
nent, appearing at the 3rd, 5th, etc. order response to 
the ijj + 2(jj field. These harmonics can be extracted by in- 
tegrating the signal over time [see Eq. I|24p ] , yielding the 
cumulative dipoles {C{t)) shown in Fig. flTC and Fig. [TT] 
C '. We denote the mean asymptotic value of (C(i)) af- 
ter the pulse as (C(oo)). As can be seen, all symmetry 
breaking effects induced by the w -f 2a; field are achieved 
while the system is being driven by the laser field. Once 
the field is turned off, the zero-frequency component dis- 
appears so that {C{t)) becomes constant and net dipoles 
no longer persist. This is because systems with discrete, 
non-degenerate, spectra cannot sustain net dipoles after 
the pulse. 

In the rigid case the zero-frequency term constitutes 
an important component of {fJ-(t)) and both fl and f4 
pulses generate a net dipole while the pulse is on. By 
contrast, when pulse f 1 is applied on the flexible chain 
(Fig. [TT] upper panels) the electron- vibrational couplings 
mute most of the effect. However, by applying a shorter 
pulse (lower panels) it is possible to generate net dipoles 
that are of the same order of magnitude as the ones ob- 
served for the rigid chain. In essence, by using shorter 
pulses one is limiting the detrimental effects that the ex- 
cited state vibronic couplings exert on the control. 

A complementary perspective on the influence of the 
vibronic couplings on the photoinduced control is ob- 
tained by considering the dependence of the symmetry 
breaking effect on the laser frequency uj. Figure [12] shows 
the asymptotic cumulative dipoles (C(oo)) observed in 



flexible and rigid chains after photoexcitation with pulses 
fl-f4 for different laser frequencies. In the rigid case, 
the uj vs. 2uj coherent control scenario is very robust, 
inducing net dipoles at most driving frequencies and 
with pulses of any time-width and intensity. The effect 
shows sharp resonances at selected frequencies. The ori- 
gin of these resonances can be identified by comparing 
the control map to the single particle spectrum of the 
rigid chain (Table HH) . Consider, for instance, the reso- 
nances displayed in Fig. [T^ at ha; —0.9, 1.18 and 1.42 
eV. At these frequencies the quantity 2uj coincides with 
the HOMO^LUMO, HOMO^ LUMO-hl (or HOMO- 
1^ LUMO), and HOMO-1^ LUMO-hl transition fre- 
quencies, respectively. By applying stronger pulses (e.g. 
Fig. [T2J3) it is possible to exploit high order multiphoton 
processes to generate dc terms, resulting in a complicated 
frequency dependence of the degree of control. 

The control map in the flexible chain is remarkably dif- 
ferent. In flexible chains, in order to generate appreciable 
symmetry breaking it is necessary to work with either suf- 
ficiently short pulses, so that only limited excited state 
dynamics occurs during the pulse (Figs. [TSD-D). or to 
apply stronger fields, as in Fig. I12B . in which the evolu- 
tion imposed by the field becomes dominant. When using 
short pulses, the resonance structure previously observed 
in the rigid-chain control case is partially maintained. 
However, in the case of long pulses the fine features ob- 
served in the rigid chain are washed away by the broad- 
ening of the electronic transitions, with only the rough 
features in the control map being maintained. Interest- 
ingly, when applying short pulses at selected frequencies 
it is possible to obtain a higher degree of symmetry break- 
ing in the flexible chain than in the rigid chain. However, 
the maximum control attainable (as a function of to) is 
always attenuated by the electron- vibrational couplings. 

Importantly, in all cases the sign and magnitude of 
symmetry breaking can be manipulated by varying the 
relative phase between the two components of the field. 
Figure [T3] exemplifies this dependence for selected fre- 
quencies and pulses. The electron- vibrational couplings 
in the system do not destroy the phase control of the elec- 
tronic dynamics, but merely attenuate its magnitude. It 
follows then that it is possible to use the techniques of 
coherent control to manipulate electronic dynamics even 
in the presence of significant vibronic couplings. The at- 
tenuation of the effect can be overcome by using stronger 
pulses or by working with laser pulses that are short or 
comparable with the time-scale for electronic dephasing. 
Similar results were found in sample studies on longer 
chains. 



IV. CONCLUSIONS 

We have investigated the possibility of inducing elec- 
tronic transport in conjugated polymers by irradiating 
the sample with u + 2uj lasers. To do so, we have fol- 
lowed the highly nonlinear dynamics of the electronic 
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and vibrational degrees of freedom in irans-polyacetylene 
oligomers in the presence oiu!+2uj laser pulses of different 
durations and intensities, and for a variety of frequencies. 
The simulations are performed in the mean field approx- 
imation, in which the nuclei are treated classically and 
the electrons quantum mechanically. The dynamics is 
followed for an ensemble of lattice initial conditions ob- 
tained by sampling the ground-state nuclear Wigner dis- 
tribution function in the harmonic approximation. 

Several important ways in which the electron- 
vibrational couplings in 7r-conjugated systems modify the 
photoinduced dynamics have been identified. The vi- 
bronic couplings introduce broadening of the electronic 
transitions and can cause ultrafast dephasing in the elec- 
tronic dynamics, as well as pronounced changes in the 
mean single-particle spectrum and IVR. 

The simulations reveal that the vibronic couplings can 
have strong detrimental effects on the control. In fact. 



it is difficult to induce and control electronic transport 
when the laser pulses are longer than the typical carbon- 
carbon vibrational period of ~ 30 fs, unless the pulse 
intensities are strong enough to completely dominate the 
dynamics. However, even in the presence of significant 
vibronic couplings, limited laser control of the electronic 
dynamics is possible via the use of pulses whose durations 
are comparable to the electronic dephasing timescale (~ 
10 fs). 

These results provide insight into the way coherent 

control works in large molecular systems with significant 
electron-vibrational couplings, and offer a characteriza- 
tion of the time-scale of decay of photoinduced breathers 
in conjugated polymers. We expect that the set of phe- 
nomena investigated here will be of generic importance 
in a wide range of materials characterized by significant 
electron-vibrational couplings. 
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